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In this series of three lectures, we discuss several aspects of high energy scattering among 
hadrons in Quantum Chromodynamics. The first lecture is devoted to a description of the 
parton model, Bjorken scaling and the scaling violations due to the evolution of parton 
distributions with the transverse resolution scale. The second lecture describes parton 
evolution at small momentum fraction x, the phenomenon of gluon saturation and the 
Color Glass Condensate (CGC). In the third lecture, we present the application of the 
CGC to the study of high energy hadronic collisions, with emphasis on nucleus-nucleus 
collisions. In particular, we provide the outline of a proof of high energy factorization 
for inclusive gluon production. 

Preprint CERN-PH-TH/2007-131 

1. Introduction 

Quantum Cliromodynamics (QCD) is very successful at describing hadronic scatter- 
ings involving very large momentum transfers. A crucial element in these successes is 
the asymptotic freedom of QCD ^ , that renders the coupling weaker as the momen- 
tum transfer scale increases, thereby making perturbation theory more and more 
accurate. The other important property of QCD when comparing key theoretical 
predictions to experimental measurements is the factorization of the short distance 
physics which can be computed reliably in perturbation theory from the long dis- 
tance strong coupling physics related to confinement. The latter are organized into 
non-perturbative parton distributions, that depend on the scales of time and trans- 
verse space at which the hadron is resolved in the process under consideration. In 
fact, QCD not only enables one to compute the perturbative hard cross-section, but 
also predicts the scale dependence of the parton distributions. A generic issue in 

* Lectures given at the Xth Hadron Physics Workshop, March 2007, Florianopolis, Brazil. 
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Fig. 1. Generic hard process in the scattering of two hadrons. Left: Leading Order. Right: Next- 
to-Leading Order correction involving gluon radiation in the initial state. 



the application of perturbative QCD to the study of hadronic scatterings is the oc- 
currence of logarithmic corrections in higher orders of the perturbative expansion. 
These logarithms can be large enough to compensate the extra coupling constant 
as they come accompanied with, thus voiding the naive, fixed order, application 
of perturbation theory. Consider for instance a generic gluon-gluon fusion process, 
as illustrated on the left of figure 1, producing a final state of momentum P^. The 
two gluons have longitudinal momentum fractions x\^2 given by 

a;i,2 = e. , (1) 



where = y P ±_ + [P = Pp.P'^ is the invariant mass of the final state) and 
Y = ln(P+/P^)/2. On the right of figure 1 is represented a radiative correction 
to this process, where a gluon is emitted from one of the incoming lines. Roughly 
speaking, such a correction is accompanied by a factor 



kl 



(2) 



where z is the momentum fraction of the gluon before the splitting, and k±^ its trans- 
verse momentum. Such corrections produce logarithms, log(l/a;i) and log(Afj^), 
that respectively become large when xi is small or when M_l is large compared 
to typical hadronic mass scales. These logarithms tell us that parton distributions 
must depend on the momentum fraction x and on a transverse resolution scale 
M_L, that are set by the process under consideration. In the linear regime'*, there 
are "factorization theorems" - fct -factorization ^ in the first case and collinear fac- 
torization ^ in the second case - that tell us that the logarithms are universal and 
can be systematically absorbed in the definition of parton distributions ^. The x de- 
pendence that results from resumming the logarithms of 1/x is taken into account 
by the BFKL equation Similarly, the dependence on the transverse resolution 
scale M± is accounted for by the DGLAP equation ^. 



''We use the denomination "linear" here to distinguish it from the saturation regime discussed 

later that is characterized by non-linear evolution equations. 

''The latter is currently more rigorously established than the former. 
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The application of QCD is a lot less straightforward for scattering at very large 
center of mass energy, and moderate momentum transfers. This kinematics in fact 
dominates the bulk of the cross-section at collider energies. A striking example of 
this kinematics is encountered in Heavy Ion Collisions (HIC), when one attempts 
to calculate the multiplicity of produced particles. There, despite the very large 
center of mass energy'^, typical momentum transfers are small*^, of the order of a 
few GeVs at most. In this kinematics, two phenomena that become dominant are 

• Gluon saturation : the linear evolution equations (DGLAP or BFKL) 

for the parton distributions implicitly assume that the parton densities in 
the hadron are small and that the only important processes are splittings. 
However, at low values of x, the gluon density may become so large that 
gluon recombinations are an important effect. 

• Multiple scatterings : processes involving more than one parton from a 
given projectile become sizeable. 

It is highly non trivial that this dominant regime of hadronic interactions is 
amenable to a controlled perturbative treatment within QCD, and the realization 
of this possibility is a major theoretical advance in the last decade. The goal of 
these three lectures is to present the framework in which such calculations can be 
carried out. 

In the first lecture, we will review key aspects of the parton model. Our re- 
curring example will bo the Deep Inelastic Scattering (DIS) process of scattering 
a high energy electron at high momentum transfers off a proton. Beginning with 
the inclusive DIS cross-section, we will arrive at the parton model (firstly in its 
most naive incarnation, and then within QCD), and subsequently at the DGLAP 
evolution equations that control the scaling violations measured experimentally. 

In the second lecture, we will address the evolution of the parton model to small 
values of the momentum fraction x and the saturation of the gluon distribution. 
After illustrating the tremendous simplification of high energy scattering in the 
eikonal limit, we will derive the BFKL equation and its non- linear extension, the 
BK equation. We then discuss how these evolution equations arise in the Color 
Glass Condensate effective theory. We conclude the lecture with a discussion of the 
close analogy between the energy dependence of scattering amplitudes in QCD and 
the temporal evolution of reaction-diffusion processes in statistical mechanics. 

The third lecture is devoted to the study of nucleus-nucleus collisions at high 
energy. Our main focus is the study of bulk particle production in these reactions 
within the CGC framework. After an exposition of the power counting rules in 
the saturated regime, we explain how to keep track of the infinite sets of diagrams 

'^At RHIC, center of mass energies range up to \/s = 200 GeV/nucleon; the LHC will collide nuclei 
at ^/s = 5.5 TeV/nucleon. 

'^For instance, in a collision at \/s = 200 GeV between gold nuclei at RHIC, 99% of the multiplicity 
comes from hadrons whose p±_ is below 2 GeV. 
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that contribute to the inclusive gluon spectrum. Specifically, we demonstrate how 
these can be resummed at leading and next-to-leading order by solving classical 
equations of motion for the gauge fields The inclusive quark spectrum is discussed 
as well. We conclude the lecture with a discussion of the inclusive gluon spectrum 
at next-to-leading order and outline a proof of high energy factorization in this 
context. Understanding this factorization may hold the key to understanding early 
thermalization in heavy ion collisions. Some recent progress in this direction is 
briefly discussed. 

2. Lecture I : Parton model, Bjorken scaling, scaling violations 

In this lecture, we will begin with the simple parton model and develop the conven- 
tional Operator Product Expansion (OPE) approach and the associated DGLAP 
evolution equations. To keep things as simple as possible, we will use Deep Inelastic 
Scattering to illustrate the ideas in this lecture. 

2.1. Kinematics of DIS 




Fig. 2. Kinematical variables in the Deep Inelastic Scattering process, k and P are known from 
the experimental setup, and k' is obtained by measuring the deflected lepton. 

The basic idea of Deep Inelastic Scattering (DIS) is to use a well understood 
lepton probe (that does not involve strong interactions) to study a hadron. The 
interaction is via the exchange of a virtual photon^. Variants of this reaction involve 
the exchange of a or boson which become increasingly important at large 
momentum transfers. The kinematics of DIS is characterized by a few Lorentz 
invariants (see figure 2 for the notations), traditionally defined as 

= P-q 

s = {P + kf 

Ml = iP + qf =ml+2v + q\ (3) 

''If the virtuality of the photon is small (in photo-production reactions for instance), the assertion 
that the photon is a "well known probe that does not involve strong interactions" is not valid 
anymore. Indeed, the photon may fluctuate, for instance, into a p meson. 
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where m^^ is the nucleon mass (assuming that the target is a proton) and the 
invariant mass of the hadronic final state. Because the exchanged photon is space- 
Uke, one usually introduces = —q^ > 0, and also x = Q'^/2u. Note that since 
M"^ > , we must have < .x < 1 - the value x = 1 being reached only in the 
case where the proton is scattered elastically. 

The simplest cross-section one can measure in a DIS experiment is the total 
inclusive electron+proton cross-section, where one sums over all possible hadronic 
final states : 

states X 

The partial cross-section associated to a given final state X can be written as 

where [d^^] denotes the invariant phase-space element for the final state X and M.^ 
is the corresponding transition amplitude. The "spin" symbol denotes an average 
over all spin polarizations of the initial state and a sum over those in the final 
state. The transition amplitude is decomposed into an electromagnetic part and a 
hadronic matrix element as 

M, = ^ [u{k')^^u{k)\ {X\J,{Q)\N{P)) . (6) 

In this equation is the hadron electromagnetic current that couples to the photon, 
and \N{P)') denotes a state containing a nucleon of momentum P. 

Squaring this amplitude and collecting all the factors, the inclusive DIS cross- 
section can be expressed as 

= -A^^L^'^W^, ' (7) 

d'^k SzTT-^^s — m^j 

where the leptonic tensor (neglecting the electron mass) is 

L"^ ^ (tl(fe')7''n(fe)«(fe)7^«(fe')),pi„ 

= 2(fc''fc'" + fc'fc'" - 5'*" fc • k') . (8) 

and W^i, - the hadronic tensor - is defined as 

AttW^, ^ J2 J[d^.]{2nfS{P + q-P,) 

states X 

X {{N{P)\4{0)\X){X\J^{0)\N{P))) 



spin 



= J d^e^'^-y {{N{P)\jt{y)MO)\N{P)))^^^^ 



(9) 



The second equality is obtained using the complete basis of hadronic states X. 
Thus, the hadronic tensor is the Fourier transform of the expectation value of the 
product of two currents in the nucleon state. An important point is that this object 
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cannot be calculated by perturbative methods. This rank-2 tensor can be expressed 
simply in terms of two independent structure functions as a consequence of 

• Conservation of the electromagnetic current : g^VF'^'' = q^W^'^ = 

• Parity and time-reversal symmetry : W^'^ = W"^ 

• Electromagnetic currents conserve parity : the Levi-Civita tensor e^^'^P'^ 
cannot appear^ in the tensorial decomposition of W^^ 

When one works out these constraints, the most general tensor one can construct 
from P'^, 9^ and g^^^ reads : 



P-q 



P-q 



Piy - qv — 5- 
1 



(10) 



where Fi^2 are the two structure functions^. As scalars, they only depend on Lorentz 
invariants, namely, the variables x and Q^. The inclusive DIS cross-section in the 
rest frame of the proton can be expressed in terms of Fi 2 as 



■N 



dE'dVL 4m„£;2sin*(0/2) 



2 Fi sin^ 



'Tit, 9 a 

—E. F2 COS^ - 

v 2 



(11) 



where represents the solid angle of the scattered electron and E' its energy. 



2.2. Experimental facts 

Two major experimental results from SLAC ^ in the late 1960's played a crucial 
role in the development of the parton model. The left plot of figure 3 shows the 




Fig. 3. SLAC results on DIS. 



measured values of F2{x,Q'^) as a function of x. Even though the data covers 
a significant range in Q^, all the data points seem to line up on a single curve, 
indicating that F2 depends very little on in this regime. This property is now 

^This property is not true in DIS reactions involving the exchange of a weak current; an additional 
structure function F3 is needed in this case. 

^The structure function F2 differs slightly from the W2 defined in : F2 = uW2/m^. 
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known as Bjorken scaling ^. In the right plot of figure 3, one sees a comparison of 
F2 with the combination^ = F2 — 2x Fi . Although there are few data points for 
F^, one can see that it is significantly lower than F2 and close to zero As we shall 
sec shortly, these two experimental facts already tell us a lot about the internal 
structure of the proton. 



2.3. Naive parton model 

In order to get a first insight into the inner structure of the proton, it is interesting to 
compare the DIS cross-section in eq. (11) and the e~/U~ cross-section (also expressed 
in the rest frame of the muon), 



da. 



e~ ix~ 



sin^ - + 



"« 2 
— cos - 



(12) 



Note that, since this reaction is elastic, the corresponding x variable is equal to 
1, hence the delta function in the prefactor. The comparison of this formula with 
eq. (11), and in particular its angular dependence, is suggestive of the proton being 

composed of point like fcrmions - named partons by Fcynman - off which the virtual 
photon scatters. If the constituent struck by the photon carries the momentum Pc, 
this comparison suggests that 



2 ~ F2 ~ (5(1 - Xc) with Xc = 



"^q-Pc 



(13) 



Assuming that this parton carries the fraction x^ of the momentum of the proton, 
i.e. Pc = XpP, the relation between the variables x and Xc is Xc = x/x^.. Therefore, 
we get : 



2Fi ~ F2 ~ XpS{x - Xp) . 



(14) 



In other words, the kincmatical variable x measured from the scattering angle of 
the electron would be equal to the fraction of momentum carried by the struck 
constituent. Note that Bjorken scaling appears quite naturally in this picture. 

Having gained intuition into what may constitute a proton, we shall now com- 
pute the hadronic tensor W'^'^ for the DIS reaction on a free fermion i carrying 
the fraction x^ of the proton momentum. Because we ignore interactions for the 
time being, this calculation (in contrast to that for a proton target) can be done in 



, the longitudinal structure function, describes the inclusive cross-section between the proton 
and a longitudinally polarized proton. 

'Prom current algebra, it was predicted that F2 = 2xFi; this relation is known as the Callan-Gross 
relation ^. 
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closed form. We obtain, 



x{{x,P\j'^\0)\p'){p'\r{0)\x,P)) 



spin 




Pf'-q' 




) 



P' 




where is the electric charge of the parton under consideration. Let us now assume 
that in a proton there are fi{Xp)dXp partons of type i with a momentum fraction 
between x^ and Xp + dx^ , and that the photon scatters incoherently off each of 
them. Wc would thus have 



(The factor x^ in the denominator is a "flux factor".) At this point, we can simply 
read the values of Fi^2, 



We thus see that the two experimental observations of i) Bjorken scaling and ii) the 
Callan-Gross relation are automatically realized in this naive picture of the proton^ . 

Despite its success, this model is quite puzzling, because it assumes that partons 
are free inside the proton - while the rather large mass of the proton suggests 
a strong binding of these constituents inside the proton. Our task for the rest 
of this lecture is to study DIS in a quantum field theory of strong interactions, 
thereby turning the naive parton model into a systematic description of hadronic 
reactions. Before we proceed further, let us describe in qualitative terms (see ^'^ for 
instance) what a proton constituted of fermionic constituents bound by interactions 
involving the exchange of gauge bosons may look like. In the left panel of figure 4 are 




Fig. 4. Cartoons of the valence partons of a proton, and their interactions and fluctuations. Left: 
proton at low energy. Right: proton at high energy. 



represented the three valence partons (quarks) of the proton. These quarks interact 




(15) 





(16) 



Jin particular, = in this model is intimately related to the spin 1/2 structure of the scattered 
partons. Scalar partons, for instance, would give F\ = 0, at variance with experimental results. 
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by gluon exchanges, and can also fluctuate into states that contain additional gluons 
(and also quark-antiquark pairs). These fluctuations can exist at any space-time 
scale smaller than the proton size 1 fcrmi). (In this picture, one should think of 
the horizontal axis as the time axis.) When one probes the proton in a scattering 
experiment, the probe (e.g. the virtual photon in DIS) is characterized by certain 
resolutions in time and in transverse coordinate. The shaded area in the picture 
is meant to represent the time resolution of the probe : any fluctuation which is 
shorter lived than this resolution cannot be seen by the probe, because it appears 
and dies out too quickly. 

In the right panel of figure 4, the same proton is represented after a boost, while 
the probe has not changed. The main difference is that all the internal time scales 
are Lorentz dilated. As a consequence, the interactions among the quarks now take 
place over times much larger than the resolution of the probe. The probe therefore 
sees only free constituents. Moreover, this time dilation allows more fluctuations 
to be resolved by the probe; thus, a high energy proton appears to contain more 
gluons than a proton at low energy"^. 



2.4. Bjorken scaling from free field theory 

We will now derive Bjorken scaling and the Callan-Gross relation from quantum 

field theory. We will consider a theory involving fermions (quarks) and bosons 
(gluons), but shall at first consider the free field theory limit by neglecting all their 
interactions. We will consider a kinematical regime in DIS that involves a large 
value of the momentum transfer and of the center of mass energy y/s of the 
collision, while the value of x is kept constant. This limit is known as the Bjorken 
limit. 

To appreciate strong interaction physics in the Bjorken limit, consider a frame 
in which the 4-momentum of the photon can be written as 

q^ = —iu,0,0jiy^ + mlQ^). (17) 

m„ V " 

From the combinations of the components of q'^ 

+ q° + q^ V 
q^ = ^— ~ > +00 

- 9° -9^ 

q = ^= — ^ ra^x constant , (18) 

v2 

and because q ■ y = q'^y~ + q~y~^ — Q± ■ y±, the integration over yf^ in W^^"^ is 
dominated by 

y"~^^0 , y+-(m„x)-i. (19) 



Equivalently, if the energy of the proton is fixed, there are more gluons at lower values of the 
momentum fraction Xp . 
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Therefore, the invariant separation between the points at which the two currents are 
evaluated is < 2y^y~ ~ — ^ 0. Noting that in eq. (9) the product of the two 
currents can be replaced by their commutator, and recalling that expectation values 
of commutators vanish for space-like separations, we also see that > 0. Thus, the 
Bjorken limit corresponds to a time-like separation between the two currents, with 
the invariant separation going to zero, as illustrated in figure 5. It is important 

Fig. 5. Region of that dominates in the Bjorken limit. 

to note that in this limit, although the invariant y^ goes to zero, the components of 
yf^ do not necessarily become small. This will have important ramifications when 
we apply the Operator Product Expansion to VK'"'. 

For our forthcoming discussion, consider the forward Compton amplitude T'^'^ 

4nT,^ ^ I J d^ye^^'-y ((iV(P)|T(jt(y) J,(0))|iV(P)))^p.^ . (20) 




Fig. 6. Forward Compton amphtude. We have also represented a cut contributing to Vy". 

It differs from W^^ by the fact that the two currents are time-ordered, and as 
illustrated in figure 6, one can recover W^^ from its imaginary part, 

W^^, = 2ImT^, . (21) 

At fixed Q^, T^'' is analytic in the variable v, except for two cuts on the real 
axis that start dX v — ±(3^/2. The cut at positive v corresponds to the threshold 
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(P + g)^ > above which the DIS reaction becomes possible, and the cut at 
negative v can be inferred from the fact that T'^" is unchanged under the exchange 
(/i y-,q ^ —q). It is also possible to decompose the tensor Ti^" in terms of two 
structure functions Ti_2 : 

V = -Ti [g,. -'-f-)+^q{P>^- 1^^) (^^ - 9^^) > (22) 

and the DIS structure functions Fi^2 can be expressed in terms of the discontinuity 
of Ti^2 across the cuts. 

We now remind the reader of some basic results about the Operator Product 
Expansion (OPE) ^^'^2. Consider a correlator (A{0)B{y)(j){xi) ■ ■ ■ (f){xn)), where A 
and B arc two local operators (possibly composite) and the cf)^s are unspecified field 
operators. In the limit 0, this object is usually singular, because products 

of operators evaluated at the same point are ill-defined. The OPE states that the 
nature of these singularities is a property of the operators A and B, and is not 
influenced by the nature and localization of the (j){xiys. This singular behavior can 
be expressed as 

A{0)B{y) = VQ(2/)0,(0), (23) 

— ^0 ' ' 

where the Ci(ij) are numbers (known as the Wilson coefEcients) that contain the 
singular y^^ dependence and the 0^(0) are local operators that have the same quan- 
tum numbers as the product AB. This expansion - known as the OPE - can then be 
used to obtain the limit y'' — » of any correlator containing the product A{0)B{y). 
If d(C'i), d(^),and d{B) are the respective mass dimensions of the operators Oi,A 
and B, a simple dimensional argument tells us that 

Ci{y) ~ |y|d(O0-dM)-d(B) (up to logarithms) . (24) 

yf^ — *0 

(Here \y\ = ^/yjly^-) From this relation, we see that the operators Oi having the 
lowest dimension lead to the most singular behavior in the limit y^ 0. Thus, only 
a small number of operators are relevant in the analysis of this limit and one can 
ignore the higher dimensional operators. 

Things are however a bit more complicated in the case of DIS, because only 
the invariant goes to zero, while the components y^ do not go to zero. The 
local operators that may appear in the OPE of T(J^(j/) Jy(0)) can be classified 
according to the representation of the Lorentz group to which they belong. Let us 
denote them Og\"''^' , where s is the "spin" of the operator (the number of Lorentz 
indices it carries), and the index i labels the various operators having the same 
Lorentz structure. The OPE can be written as : 



(25) 
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Because they depend only on the 4- vector y^^, the Wilson coefScients must be of 
the form' 

ct;:■■■^M^v,.■■■y^sCs.^{y'). (26) 

where Cs^iiy"^) depends only on the invariant y^. Similarly, the expectation value of 
the operators O'^'^'"'^' in the proton state can only depend on the proton momentum 
P'', and the leading part in the Bjorken limit is™ 

((iV(P)|Of/-^'(0)|7V(P))\ =P'^i---P''^ , (27) 

\ ' / spin 

where the (Os,i) are some non-perturbative matrix elements. 

Let us now denote by ds,i the mass dimension of the operator O^y Then, 
the dimension of Cs,i{y'^) is 6 + s — ds,i, which means that it scales like 

CsM) ~ (j;2)(d,.-«-6)/2 _ (28) 

Because the individual components of y^^ do not go to zero, it is this scaling alone 
that determines the behavior of the hadronic tensor in the Bjorken limit. Contrary 
to the standard OPE, the scaling depends on the difference between the dimension 
of the operator and its spin, called its twist ts^i = dgs — s, rather than its dimension 
alone. The Bjorken limit of DIS is dominated by the operators that have the lowest 
possible twist. As we shall see, there is an infinity of these lowest twist operators, 
because the dimension can be compensated by the spin of the operator. If we go 
back to the structure functions Ti_2, we can write 

Tr{x, Q2) = ^ x'^^-' (^«.^) I^r;sAQ^) {v = 1,2), (29) 

s i 

where ai = and 02 = 1. The difference by one power of x (at fixed Q^) between 
Ti and T2 comes from their respective definitions from T'^" that differ by one power 
of the proton momentum P. Eq. (29) gives the structure functions Ti^2 as a series 
of terms, each of which has factorizcd x and dependences. (The functions Dr-^s.i 
(r = 1,2) are related to the Fourier transform of Cs,i{y'^), and thus can only depend 
on the invariant Q'^). Moreover, for dimensional reasons, the functions Dr-s,i must 
scale like Q^^'*"'*'-'. Therefore, it follows that Bjorken scaling arises from twist 2 
operators. It is important to keep in mind that in cq. (29), the functions Dr-s,i are 
in principle calculable in perturbation theory and do not depend on the nature of 
the target, while the (Os,i)'s are non perturbative matrix elements that depend on 
the target. Thus, the OPE approach in our present implementation cannot provide 
quantitative results beyond simple scaling laws. 

'There could also be terms where one or more pairs yi^'y^i are replaced by g^*^i , but such 
terms axe less singulars in the Bjorken limit. 

™Here also, there could be terms where a pair pMipMj jg replaced by g^*^^ , but they too 
lead to subleading contributions in the Bjorken limit. 
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It is easy to check that Ti is even in x while T2 is odd; this means that only 
even values of the spin s can appear in the sum in eq. (29). We shall now rewrite 
this equation in a more compact form to see what it tells us about the structure 
functions Fi.2- Writing 

Tr^ tris.Q^) X'^"-^ = tr{s,Q') (-^) " l^'-'^^ , (30) 

even s even s \ ^ / 

we get (for 5 even) 

where C is a small circle around the origin in the complex v plane (see figure 7). This 




Fig. 7. Contour in the complex i/ plane, and its deformation to pick up the contribution of the 
cuts. 



contour can then be deformed and wrapped around the cuts along the real axis, as 
illustrated in the figure 7. Because the structure function Fr is the discontinuity of 
Tr across the cut, we can write 

U{s, Q2) ^1 f'^ a;--. Fr{x, Q2) . (32) 

Therefore, we see that the OPE gives the x-moments of the DIS structure functions. 

In order to go further and calculate the perturbative Wilson coefficients Dr;s,i, 
we must now identify the twist 2 operators that may contribute to DIS. In a theory 
of fermions and gauge bosons, we can construct two kinds of twist 2 operators : 

Qi^r-t^. = Fj'^'d''^ ■ . , (33) 

where the brackets {• • • } denote a symmetrization of the indices /ii • • • /i^ and a 
subtraction of the traced terms on those indices. To compute the Wilson coefficients, 
the simplest method is to exploit the fact that they are independent of the target. 
Therefore, we can take as the "target" an elementary object, like a quark or a 
gluon, for which everything can be computed in closed form (including the (O^.i)). 
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Consider first a quark state as the target, of a given flavor / and spin a. At lowest 
order, one has 

(/,a|O,^--^=|/,^)=0. (34) 

Averaging over the spin, and comparing with P^^ • • • P'^' {Os,i), we get 

(a,/')^ = <5//' , {Os,a)j = ^. (35) 

On the other hand, we have already calculated directly the hadronic tensor for a 
single quark. By computing the moments of the corresponding Fi^2, we get the 
tr{s, Q"^) for s even : 

<i(s,Q2)^ie2 ^ ^{s,Q^) = -e'f . (36) 
Prom this, the bare Wilson coefficients for the operators involving quarks are 

i?w(Q') = -e^ , D2;sj{Q^) = -e) . (37) 

By repeating the same steps with a vector boson state, those involving only gluons 
are 

Dx.sAQ'') = D2;s,g{Q^) = , (38) 

if the vector bosons are assumed to be electrically neutral. 

Going back to a nucleon target, we cannot compute the {^Os.i}- However, we can 
hide momentarily our ignorance by defining functions // (a;) and ff{x) (respectively 
the quark and antiquark distributions) such that" 



I 



^ dx , 



ff{x)+ff{x) ={Osj) . (39) 



(The sum ff{x) + ff{x) is known as the singlet quaj-k distribution ofBavorf.) Thus, 
the OPE formulas for Pi and P2 on a nucleon in terms of these quark distributions 
are 

^i(^) = jE4f//(^)+//M ' F2{x)=2xF,{x). (40) 

We see that these formulas have the required properties: (i) Bjorken scaling and 
(ii) the Callan-Gross relation. 

Despite the fact that the OPE in a free theory of quarks and gluons leads to a 
result which is embarrassingly similar to the much simpler calculation we performed 
in the naive parton model, this exercise has taught us several important things : 



°DIS with exchange of a photon cannot disentangle the quarks from the antiquaries. In order to 
do that, one could scatter a neutrino off the target, so that the interaction proceeds via a weak 
charged current. 
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• We can derive an operator definition of tlie parton distributions fi{x) (al- 
beit it is not calculable perturbatively) 

• Bjorken scaling can be derived from first principles in a field theory of 
free quarks and gluons. This was a puzzle prc-QCD because clearly these 
partons are constituents of a strongly bound state. 

• The puzzle could be resolved if the field theory of strong interactions be- 
came a free theory in the limit Q"^ +00, a property known as asymptotic 
freedom. 

As shown by Gross, Politzer and Wilczek in 1973, non-Abclian gauge theories with 
a reasonable number of fermionic fields (e.g. QCD with 6 flavors of quarks) are 
asymptotically free^ and were therefore a natural candidate for being the right 
theory of the strong interactions. 

2.5. Scaling violations 

Although it was interesting to see that a free quantum field theory reproduces the 
Bjorken scaling, this fact alone does not tell much about the detailed nature of the 
strong interactions at the level of quarks and gluons. Much more interesting are 
the violations of this scaling that arise from these interactions and it is the detailed 
comparison of these to experiments that played a crucial role in establishing QCD 
as the theory of the strong interactions. 

The effect of interactions can be evaluated perturbatively in the framework of 
the OPE, thanks to renormalization group equations. In the previous discussion, 
we implicitly assumed that there is no scale dependence in the moments (Cs,i) of 
the quark distribution functions. But this is not entirely true; when interactions 
are taken into account, they depend on a renormalization scale jj^. The parton 
distributions become scale dependent as well. However, since Fi^2 are observable 
quantities that can be extracted from a cross-section, they cannot depend on any 
renormalization scale. Thus, there must also be a dependence in the Wilson co- 
efficients, that exactly compensates the /i^ dependence originating from the (Os,i)- 
By dimensional analysis, the Wilson coefficients have an overall power of set by 
their dimension (sec the discussion following eq. (29)), multiplied by a dimension- 
less function that can only depend on the ratio Q"^ j \j? . By comparing the Callan- 
Symanzik equations^^ for T'''' with those for the expectation values {Os,i), the 
renormalization group equation^^ obeyed by the Wilson coefficients is ° 

[{-Qd^ + (3{g)dg) 6ij - 7.,,i(5)] £>r;.,i (Q/m, 9) = , (41) 

where (3{g) is the beta function, and ^ajiig) is the matrix of anomalous dimensions 
for the operators of spin s (it is not diagonal because operators with identical 
quantum numbers can mix through renormalization). 

°We have used the fact that the electromagnetic currents are conserved and therefore have a 
vanishing anomalous dimension. Note also that we have exploited the fact that for twist 2 operators 
Dj.-sj depends only on Q'^ j ij? , so that we can replace pQ^j, by —QBq. 
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In order to solve these equations, let us first introduce the running coupling 
g{Q,g) such that 

'9 



HQ/Qo) = f 

J a 



(42) 



Note that this is equivalent to QdQg{Q,g) = (3{g{Q,g)) and g{Qo,g) = g', in other 
words, g{Q, g) is the value at the scale Q of the coupling whose value at the scale Qo 
is g. The usefulness of the running coupling stems from the fact that any function 
that depends on Q and g only through the combination g{Q,g) obeys the equation 



[-QdQ+mdg]F{g(Q,g))^0 



(43) 



It is convenient to express the Wilson coefficients at the scale Q from those at the 
scale Qo as 



Dr-sAQIiJ'^g) = Dr;s,j{Qo/iJ',g{Q,g)) 



e-/Qo'^7»(s(M,s)) 



(44) 



In QCD, which is asymptotically free, we can approximate the anomalous dimen- 
sions and running coupling at one loop by 



is,ijig) = g Aij{s) , g {Q,g) 



(45) 



(The Aij (s) are obtained from a 1-loop perturbative calculation.) In this case, the 
scale dependence of the Wilson coefficients can be expressed in closed form as 



-Dr;s,i(<9//^,5) = Dr;s,jiQo/iJ',giQ,g)) 



( HQ/ 



Vln(go/AQCD) 



^0 



A{s) 



(46) 



From this formula, we can write the moments of the structure functions. 



f HQ/K 



{Os 



Qo 



(47) 



J fi 



(and a similar formula for F2). Wc sec that we can preserve the relationship between 
Fi and the quark distributions, eq. (40), provided that we let the quark distributions 
become scale dependent in such a way that their moments read 



dx 



/ ln(Q/A,,,J 
\ln{Qo/AQcn) 



'-A(s 



iOs,)^^ .(48) 

fi 

By also calculating the scale dependence of F2, one could verify that the Callan- 
Gross relation F2{x, Q"^) = 2xFi{x, Q^) is preserved at the 1-loop order. It is crucial 
to note that, although we do not know how to compute the expectation values 
{Os,i) at the starting scale Qo, QCD predicts how the quark distribution varies 
when one changes the scale Q. We also see that, in addition to a dependence on Q^, 
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the singlet quark distribution now depends on the expectation value of operators 
that involve only gluons (when the index i = g in the previous formula). 

The scale dependence of the parton distributions can also bo reformulated in 
the more familiar form of the DGLAP equations. In order to do this, one should 
also introduce a gluon distribution fg, also defined by its moments, 



Jo 



/ ln(Q/A, 



\\n{Qo/AQCD) 



f>0 



A{s) 



{Os,i) 



Qo 



(49) 



Then one can check that the derivatives of the moments of the parton distributions 
with respect to the scale are given by 



Q 



Aji{s)fj{s,Q^) , (50) 
In order to turn this 



where we have used the shorthands // = // + fj , fg = fg 
equation into an equation for the parton distributions themselves, one can use 

'■^ dx „ dy 



A{s)f{s)= [ -x^ [ ^Aix/y)f{y) 
Jo X y 



(51) 



that relates the product of the moments of two functions to the moment of a 
particular convolution of these functions. Using this result, and defining splitting 
function Pij from their moments. 



X ^ij (*^) 

X 



-4n'^Aij{s) , 



it is easy to derive the DGLAP equation^. 



Q 



^dfi{x,Q^) _gHQ,9) f'dy 



87r2 



■Mx/y)fJy,Q^ 



(52) 



(53) 



that resums powers of aslog{Q'^ /Qq). This equation for the parton distributions 
has a probabilistic interpretation : the splitting function g'^Pji{z) \n{Q-^) can be seen 
as the probability that a parton j splits into two partons separated by at least 
(so that a process with a transverse scale Q will see two partons) , one of them being 
a parton i that carries the fraction z of the momentum of the original parton. 
At 1-loop, the coefficients Aij{s) in the anomalous dimensions are 



1 

12' 



1 



s(s-l) (s+l)(s+2) 
2 1 



" 1 




s(s + l)(s + 2) 
2 




(54) 
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where Nf is the number of flavors of quarks. On can note that, since Agf{s) is fla- 
vor independent, the non-singletP linear combinations CLf^sj with J2f'^f — 0) 
are eigenvectors of the matrix of anomalous dimensions, with an eigenvalue A//(s). 
These linear combinations do not mix with the remaining two operators. ^ s.f 

and Cs,g, through renormalization. By examining these anomalous dimensions for 
s = 1, we can see that the eigenvalue for the non-singlet quark operators is vanish- 
ing : A^f[s = 1) = 0. Going back to the eq. (50), this implies that 



d 



/ 



(55) 



for any linear combination such that ay = 0. This relation implies for instance 
that the number of m + u quarks minus the number of d + d quarks does not depend 
on the scale Q, which is due to the fact that the splittings g ^ qq produce quarks 
of all flavors in equal numbers (if one neglects the quark masses). An interesting 
relation can also be obtained for s = 2. For this moment, the matrix of anomalous 
dimensions in the singlet sector. 



Aff{2) Afg{2}\ _ 

NfAgji2) Agg{2) 



1 




(56) 



has a vanishing eigenvalue, which means that a hnear combination of the flavor 

•J- 



singlet operators is not renormalized : Oj.g + J2f ^2j- ^^^^ leads also to a sum 



rule 



dxx 



[ff{x,Q') + ff{x,Q^)\+fg{x,Q^) 

f 



(57) 



whose physical interpretation is the conservation of the total momentum of the 
proton - which therefore cannot depend on the resolution scale Q. (CoUinear split- 
tings, that are responsible for the Q dependence of the number of partons, do not 
alter their total momentum.) 

We have seen that QCD can be used to calculate the value of the Wilson coeffi- 
cients as well as the scale dependence of the non-perturbative parton distributions. 
In practice, when one compares DIS data with theoretical predictions, one needs 
only to adjust the value of the parton distributions at a relatively low initial scale 
Qo, and then one uses the DGLAP evolution equations in order to obtain their value 
at a higher Q. This program has now been implemented to three loops (NNLO), 
and has been very successful in explaining the inclusive DIS data. The agreement 
between QCD and the DIS measurements is illustrated in figure 8 (see for instance 
for more details). 



PHere, the word "singlet" refers to the flavor of the quarks. 
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Fig. 8. Comparison of the measured F2 with QCD fits. 



3. Lecture II : Parton evolution at small x and gluon saturation 

In the first lecture, we introduced the parton model and the evolution of parton 
distributions with the transverse resolution scale - and the corresponding resum- 
mation of the powers of as log(Q^). We now turn to the logarithms of l/x. These 
logarithms are expected to be the dominant effect in processes where the collision 
energy ^/s is much larger than the typical transverse momentum scale involved in 
the process, and may lead to gluon saturation at very small x. 



3.1. Eikonal scattering 

Before going to the main subject of this lecture, let us make a detour through an 
important result concerning the high energy limit of the scattering amplitude of 
some state off an external field. Our derivation here follows Consider the generic 
S'-matrix element 



Sf3a = (/3out|ain) = ( An | C^(+00 , -Oo) | ttin) 

for the transition from a state a to a state /3 where 



U{+oo,—(X)) = r+ exp i d a; £int(0in(a;)) 



(58) 



(59) 



is the evolution operator from t = —00 to t = +00. (r+ denotes an ordering in 
the light-cone time .) The interaction Lagrangian Ci^t contains both the self- 
interactions of the fields and their interactions with the external field. Now apply 
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a boost in the z direction to all the particles contained in the states a and /3. 
Formally, this can be done by multiplying the states by exp(— iwif^), where a; is 
the rapidity of the boost and the generator of longitudinal boosts. Our goal is 
to compute the limit u) — > +00 of the transition amplitude, 

SfJ= lim (A„|e^-^'[/(+oo,-oo)e-^-^>i„) . (60) 

^ W— >+(X) ^ ' ' ' 

The behavior of scattering amplitudes in this limit is easy to understand. The time 
spent by the incoming particles in the region where the external field is acting goes 
to zero as the inverse of the collision energy E. If the coupling to the external field 
was purely scalar, this would imply that the scattering amplitude itself goes to zero 
as . However, in the case of a vector coupling, the longitudinal component of 
the current increases as i?, which compensates the decrease in the interaction time, 
thereby leading to a finite (non-zero and non infinite) high energy limit. 

For this reason, let us assume that the coupling of the fields to the external 
potential is of the form g Ay^{x)J^[x) where is a vector current built from the 
elementary fields of the theory under consideration. In order to simplify the discus- 
sion, we also assume that the external potential is non-zero only in a finite range 
in a;+, a;+ e [— i, -\-V\ (this is to avoid complications with long range interactions). 
The action of on states and operators is 

e--^'aL((z)e-^' = aL(e-g+, e^g", q^) 
e-^'^^' |p- • -in) = |(e>+,p^) •••;„) 

e''^^Vin(x)e-^'^^' = </)i„(e-'^x+,e'^x-,a;i) , (61) 

namely, it multiplies the + component of momenta by e*^ and their minus com- 
ponent by e~", while keeping the transverse components unchanged. The external 
potential A^{x) is unaffected by K"^, and the components of J^{x) are changed as 
follows: 

e"'^ J-(e-"x+,e'^x-,Xi) 
j+(e-'^x+,e'"a;-,a;_L) 

Because does not modify the ordering in a;+, we can write 

e''^-^'c/(+oo, -oo)e-^"^' = T+ expi ^ d^x Ant(e'"^Vin(x)e-*'"^') . (62) 

In addition, we can split the evolution operator into three factors 

[/(-Foo, -00) = ?7(-Foo, -FL)f/(-FL, -L)U{-L, -00) (63) 

so that only the factor in the middle contains the external field. In order to deal 
with the first and last factor after the boost, it is sufficient to change variables 



e*'^^'j-(x)e'*'^^' = 
e*"^'j+(x)e-«^^' = 



21 



lim e 
lim e 



a; . This leads to 



iuiK- 



U{+oo,+L)e~ 
U{—L, — oo)e' 



f/o(+oo,0) 
= C/o(0,-oo) , 



(64) 



where Uq is the same as U, but with the self-interactions only. For the factor 
U{L, —L), the change of variables e'^x^ x^ gives us 



lim e 



iojK' 



U{+L,-L)e 



-iujK-' 



= T+ exp 



igj 



with 



x{x±) = j dx^ A (a;+,0,x_L), 



(65) 



(66) 



dx J^(0, a; , x±_) . 



Only the minus component of the external vector potential matters, because this 
is the component that couples to the longitudinal current J"*" which is enhanced 
by the boost. Therefore, the high energy limit of the transition amplitude can be 
written as 



q(00) 

'^0a 



(An|t/o(+oo,0) T+ exp ig j x{xi^)p{xa_) ;7o(0, -00) |ai„) 



(67) 



This limit is known as the eikonal limit. It is important to keep in mind that this 
formula is the exact answer for the high-energy limit; no perturbative expansion 

has been made yet, and the formula still contains the self- interactions of the fields 
of the theory to all orders. A remarkable feature of eq. (67) is that it separates the 
self-interactions of the fields and their interactions with the external potential in 
three different factors, a property which is strongly suggestive of the factorization 
between the long and short distance physics in high energy hadronic interactions. 

In order to use eq. (67) in practice, it is necessary to make an expansion in the 
self-interactions of the fields, by introducing complete sets of states between the 
three factors, 

x(7in|T+ exp ig j x{xi_)p{xi_) |(5in)((5in|J7o(0, -oo)|ai„) . (68) 



The factor |5i„)((5i„|C/(0, — oo)|ain) is the Fock expansion of the initial state. 
It reflects the fact that the state a prepared at a;+ = —00 may have fluctuated 
into another state S before it interacts with the external potential. There is also a 
similar expansion for the final state. Assuming that we have performed the Fock 
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expansion to the desired order'!, one needs to evaluate matrix elements such as 

(7in|exp ig j Xa{xi_)p''{xi_) |(5i„) . (69) 

We have reinstated color indices in this formula, since we have applications to QCD 
in mind. In order to calculate this matrix clement, the first step is to express the 
operator p°'{x±_) in terms of creation and annihilation operators of the particles 
that can couple to the external potential. For instance, the contribution that comes 
from the quarks and antiquarks is given by 

-4(p+,p±;*)rfin(p+,'7±;j>-'<^--''-^'""} • (70) 

(The quarks come with a positive sign and the antiquarks with a negative sign.) 
The contribution of the gluons would be similar, but the color matrix would be 
replaced by an clement of the adjoint representation. From this formula, we see 
that in eq. (69), the states 5 and 7 must have the same particle content, because 
each annihilation operator in is immediately followed by a creation operator 
that creates a particle of the same nature. The + component of the momenta of 
the particles in 5 and 7 must also be identical. The only difference between the 
states 5 and 7 is in the transverse momenta and in the color of their particles. 
In order to recover the eikonal limit in a more familiar form, one should go to 
impact parameter representation by performing a Fourier transformation of all the 
transverse momenta in the intermediate states 5 and 7, by defining the light-cone 
wavefunction 

^Sa{{k+,Xi^}) ^ TT / ^e-^'^-— <(5in|C/o(0,-Oo)|ain> . (71) 

Then, from the explicit form of p", it is easy to check that the only effect of the 
external potential is to multiply the function by a phase factor for each particle 
in the intermediate state : 

^Sa{{kt,Xi±}) > '^Sa{{kt,Xi±})'[[Ui{x_L) 

UiXxjL) ^T+exp ig^ J dx+ Aaix+,0,x^)t'' . (72) 

In the case of non-abelian interactions, these phase factors Ui{x±) are known as 
Wilson lines. Wilson lines resum multiple scatterings off the external field, as one 
can see by expanding the exponential. Thus, the physical picture of high energy 
scattering off some external field is that the initial state evolves from —00 to 0, 

^The main difference compared to the usual perturbation theory is that the integrations over a;+ 
run only over half of the real axis, e.g. [— oo,0]. In Fourier space, this implies that the minus com- 
ponent of the momentum is not conserved at the vertices, and that one gets energy denominators 
instead of delta functions. 
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multiply scatters during an infinitesimally short time off the external potential, 
and evolves again from to +00 to form the final state, as illustrated in figure 9. 
In terms of light-cone wavefunctions and of Wilson lines, the high energy limit of 




Fig. 9. Scattering off an external potential in the high energy limit. 



the transition amplitude reads 
dk+ 



^(oo) 
'/3q 



E 

s 



[n 



:d Xi± 



3.2. BFKL equation 

Let us now derive the BFKL equation. Our derivation is inspired from ^'^ . 
Consider the forward scattering off an external field of a state a whose simplest Fock 
component is a color singlet quark-antiquark pair. Thus, the transition amplitude 
can be written as 



i^[U{x^)U\y^)\ . (74) 



We will not need to specify more the light-cone wavefunction of the state under 
consideration. Note that the product of the two Wilson lines is traced, because the 
state a is color singlet. A crucial property of this transition amplitude is that it 
is completely independent of the collision energy. However, as we shall see, a non 
trivial energy dependence arises in this amplitude because of large logarithms in 
loop corrections. 

Consider now the 1-loop corrections to this amplitude depicted in figure 10. 
These 1-loop corrections all involve one additional gluon attached to the quark or 
antiquark lines. In some of the corrections, that we shall call real corrections^ the 
gluon is present in the state that goes through the external field. In the other cor- 
rections, the virtual corrections^ the gluon is just a fluctuation in the wavefunction 
of the initial or final state. The calculation of these diagrams is straightforward 
in the impact parameter representation. One simply needs the formula for the qqg 
vertex : 



^f(")( 



(75) 
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<$> <$> <f> <$> 

<|> <|> <z|> - - 



Fig. 10. One-loop corrections to the scattering of a dipole off an external field. Only half of the 
virtual corrections have been represented. 



where is the polarization vector of the gluon and k± its transverse momentum, 
and its expression in impact parameter space, 



ik±- (x±—z 



Armed with these tools, it is easy to obtain expressions such as 



(76) 



(fz± {x± - z±) ■ {x± - z±) 



2a, 



and 



<zr> 



x4«« / — 



ti-[t-uix^)uHyi_)t''] 

d^z± {x± - z±) ■ {yj_ - z±) 
(27r)2 {x^^z^)^{y^-z^y 

We find that the sum of all the virtual corrections reads 



Cftts 



dk+ 



(77) 



(78) 



{x±_ - zj_)'^{yj_ - zj_) 



¥'>\x^,y^) tr[U{xi_)U\y^)] , 

(79) 



where Cf = t^'t'^ = {N^ ~ 1)/2N for SU(N). In this formula, k+ is the longitu- 
dinal momentum of the gluon. As one can see, there is a logarithmic divergence 
in the integration over this variable. The lower bound should arguably be some 
non-perturbative hadronic scale A, and the upper bound must be the longitudinal 
momentum of the quark or antiquark that emitted the photon. Hence we have 
a log(p+/A), which is a large factor in the limit of high-energy (strictly speaking, 
the high-energy limit is ill defined because of these corrections) . The calculation of 
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the real corrections is a bit more involved. For instance, one has 

f dk+ fd^zj^- {xj_ - zj_) ■ {xj_- zj_) 

"^"^ J— J wr'^^'^'^ hx^-z^nx^-z^r ' ^'"^ 

where Uabi^i) is a Wilson line in the adjoint representation that represents the 
eikonal phase factor associated to the gluon {z± is the impact parameter of the 
gluon). In order to simplify the real terms, we need the following relation between 
fundamental and adjoint Wilson lines, 

eUab{Zi_)^U{Zi_)t''U\z^) , (81) 

and the Fierz identity obeyed by fundamental SU(N) matrices : 

t\Ai ^ - ^S^jSkl . (82) 

Thanks to these identities, one can rewrite all the real corrections in terms of 
the quantity S{x±,yj_) = tr [U{x±)W {yj_)'\ /N . Collecting all the terms, and 
summing real and virtual contributions, we obtain the following expression for the 
1-loop transition amplitude 

I J {xj_ - zj_y{yj_- z±_y 



27r2 



:[s{x^,y^) - S{x^,z^)S{z^,y^)} , (83) 



where we denote Y = ln(p+/A). This correction to the transition amplitude is not 
small when aj^ < Y, which means that n-loop contributions should be considered 
in order to resum all the powers (a^l")". Here, we are just going to admit that this 
n-loop calculation amounts to exponentiating the 1-loop result. In other words, 
eq. (83) is sufficient in order to obtain the derivative dS/dY, 

dS{xj_,yj_) a^Nc f ^ {x±-y±f 

a z± 



dY 2^2 J - [x^-z^Y{y^-z^Y 

x[s{x^,y^) - S{x^,z^)S{z^,y^)] . (84) 

It is customary to rewrite this equation in terms of T- matrix elements, T{x±^ yj_) = 
1 — S{x±,yj_). The BFKL equation'* describes the regime where T{x±,yj_) is small, 
so that we can neglect the terms that are quadratic in T. It reads : 

dT{x^,y^) _a,N, f ^^^^ {x^-y^f 



dY 2^2 y - — {x^-z^y{y^-zi_y 

x{T{x^,z^)+T{z^,y^^)-T{x^,y^^)} . (85) 

One can verify easily that T = is a Exed point of this equation (the right hand 
side vanishes if one sets T ~ Q), but that this fixed point is unstable (if one sets 
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T = e > 0, the right hand side is positive). Since there are no other fixed points, 
solutions of the BFKL have an unbounded growth in the high energy Umit {Y — > 
+00). This behavior however is not physical, because the unitarity of scattering 
amplitude implies that T{x±,y±) should not become greater than unity. 

3.3. Balitsky-Kovchegov equation 

The solution to the above problem was in fact already contained in eq. (84). When 
written in terms of T without assuming that T is small, 

x{t{x±,z±) +T{z±,y^) -T{x±,y^) -T{x±,z±)T{z±,y^)} , (86) 

it has a non-linear term that confines T to the range [0, 1]. Indeed, the presence of 
this quadratic term makes T = 1 a stable fixed point of the equation. Therefore, 
the generic behavior of solutions of cq. (86) is that T starts at small values at small 
Y and asymptotically reaches the value T = 1 in the high energy limit. Eq. (86) is 
known as the Balitsky-Kovchcgov cquation^"^'^^ . 

The interaction of a color singlet dipole with an external color field is a possible 
description of DIS, in a frame in which the virtual photon splits into a quark- 
antiquark pair long before it collides with the proton (the external color field would 
represent the proton target). Although it is legitimate to treat the proton as a 
frozen configuration of color field due to the brevity of the interaction, we do not 
know what this field is. Moreover, since this field is created by the partons inside the 
proton, that have a complicated dynamics, this color field must be different for each 
collision, and should therefore be treated as random. Therefore, in order to turn our 
dipole scattering amplitude into an object that we could use to compute the DIS 
cross-section at high-energy, we must average over all the possible configurations of 
the external field. Let us denote by ( • • • ) this average. The effect of this average on 
the energy dependence of the amplitude is simply taken into account by taking the 
average of eq. (86). However, one sees that the evolution equation for (T'^ involves 
in its right hand side the average of a product of two T's, (TT). Therefore, we 
do not have a closed equation anymore. An evolution equation for (TT) could be 
obtained by the same procedure, which would depend on yet another new object, 
and so on. At the end of the day, one in fact obtains an infinite hierarchy of nested 
equations, known as Balitsky's equations^^. 

It is only if one assumes that the averages of products of amplitudes factorize 
into products of averages, 

(TT) « (T) (T) , (87) 

that this hierarchy can be truncated into a closed equation which is identical to 
eq. (86) - the BK equation - with T replaced by (t). This approximation amounts 
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to drop certain correlations among the target fields, and is believed to be a good 
approximation for a large nucleus in the limit of a large number of colors^^. 

3.4. Gluon saturation and Color Glass Condensate 

The problem encountered with the indefinite growth of the solutions of the BFKL 
can be understood in terms of the behavior of the gluon distribution at small 
momentum fraction x. Indeed, in the regime where the dipole scattering amplitude 
T is still small, it can be calculated perturbatively, 

T{x^,y^) (X \x^-y^\^ xG{x,\x^ - y^\-^) , (88) 

where Y = ln(l/a;). This formula is an example of the duality that exists in the 
description of scattering processes at high energy. In the derivation of the BFKL 
and BK equations, we have treated the proton target as given once for all, and the 
energy dependence has been obtained by applying a boost to the dipole projectile. 
But, thanks to the fact that transition amplitudes are Lorentz invariant quantities, 
they can also be evaluated in a frame where the dipole is fixed, and the boost 
is applied to the proton. In this frame, the energy dependence of the scattering 
amplitude comes from the x dependence of the proton gluon distribution. 

Thus, an exponential behavior of T is equivalent to an increase of the gluon 
distribution as a power of 1/x : 

T-e'^^ < — > xG{x,Q^)^^. (89) 

(This growth of the gluon distribution is due to gluon splittings.) However, the 
gluon distribution cannot grow at this pace indefinitely. Indeed, at some point, the 
occupation number of the gluons will become large and the recombination of two 
gluons - not included in the BFKL equation - will be favored. This phenomenon is 
known as gluon saturation^'^ . In the linear regime, described by the BFKL equation. 




Fig. 11. Gluon saturation : merging of the gluons ladders initiated by two valence partons. The 
proton target is at the top of the picture and the probe at the bottom. 



each valence parton from the proton initiates its own gluon ladder (see figure 11) 
that evolves independently from the others. In the saturated regime, these gluon 
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ladders can merge, thereby reducing the growth of the gluon distribution. The effect 
of these recombinations on the scattering ampUtude is taken into account by the 
non-hnear term of the BK equation. 

A semi quantitative criterion for gluon saturation can be obtained^" by com- 
paring the surface density of gluons, p ~ xG{x, Q^)/TrR'^, and the cross-section for 
gluon recombination, a ~ ag/Q^- Saturation occurs when 1 < pa, i.e. when 

Q^<Ql , withQ? ^ a.^GMl) ^ ^r/a A_ . (90) 

A 

The quantity Qg is known as the saturation momentum. Its dependence on the 
number of nucleons A (in the case of a nuclear target) comes from the fact that 
xG{x, Q^) scales like the volume, while ttR^ is an area. Its x dependence is a phe- 
nomenological parameterization inspired by from fits of HERA data. From eq. (90), 




one can divide the x, in two regions, as illustrated in figure 12. The saturated 
regime corresponds to the domain of low Q and low x. 

Although the BK equation describes the evolution of the dipole scattering ampli- 
tude into the saturation regime, there is an equivalent description of this evolution 
- the Color Glass Condensate - in which the central role is played by the target. 
The CGC description divides the degrees of freedom in the proton into fast partons 
(large x) and slow partons (small x)^^. The fast partons are affected by time di- 
lation, and do not have any significant time evolution during the brief duration of 
the collision; therefore, they are treated as static objects that carry a color source. 
These color sources produce a current, 

jt^ = S''+S{x-)p{xi_) , (91) 

written here for a projectile moving in the +z direction. The function p{xj_) de- 
scribes the distribution of color charge as a function of the impact parameter. 
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The slow partons, on the other hand, have a non trivial dynamics during the col- 
lision, and must be treated as gauge fields. The only coupling between the fast 
and slow partons is a coupling Af^J^ between the color current of the fast partons 
and the gauge fields, which allows the fast partons to radiate slower partons by 
bremsstrahlung. Because the configuration of the fast partons prior to the collision 
is different in every collision, the function p{x±) must be a stochastic quantity, 
for which one can only specify a distribution [p] ■ Observables like cross-sections 
must be averaged over all the possible configurations of p with this distribution. In 
fact, in the CGC description, this averaging procedure is equivalent to the target 
average of the scattering amplitude that was introduced in the discussion of the 
BK equation, 

(•••)^ I [Dp]WAp] ■■■ . (92) 

A crucial point is that the distribution ^^[p] depends on Y, the rapidity that 
separates what is considered fast and slow. Because such a separation is arbitrary, 
physical quantities cannot depend on it; one can derive from this requirement a 
renormalization group equation for Wy [p] - known as the JIMWLK equation^"^ -, 
of the form : 

^^ = n[p]wAp]. (93) 

The JIMWLK Hamiltonian Hip] contains first and second derivatives with respect 
to the source p, 

n[p]= [ (^{x±)Y-r^ + l [ Xix±,y^) ^ , (94) 

where a{x±) and x(a;j^,yj_) are known functionals of p. In fact, the JIMWLK 
equation is equivalent to the infinite hierarchy of Balitsky's equations - of which the 
BK is an approximation that neglects some correlations. In the CGC description 
of scattering processes, the energy dependence of amplitudes arises from the Y 
dependence of the distribution Wy. [p] . For instance, the dipole scattering amplitude 
would be written as 



(r(x^,y^))= J [Dp] Wy[p] l-l-tT{U{x^)U^{y^)) 



(95) 



where the Wilson line U is evaluated in the color field generated by the configu- 
ration p of the color sources. This formula is very similar - at least in spirit - to 
the standard coUinear factorization in DIS. The functional [p] can be seen as 
an extension of the usual concept of parton distribution, that contains information 
about parton correlations beyond the mere number of partons, while the square 
bracket is the analogue of the "perturbative cross-section" . This formula is a Lead- 
ing Logarithm (LL) factorization formula in the sense that it resums all the powers 
(a^y)". Moreover, it also resums all the rescattering corrections, in {Qs/pi^Y, a 
feature which is not included in coUinear factorization. 
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Eq. (93) predicts the energy dependence of the distribution of sources. However, 
it must be supplemented by an initial condition at some Yq. As with the DGLAP 
equation, the initial condition is non-perturbative, and one must in general model 

it or guess it from experimental data. In the case of large nuclei, one often uses the 
McLerran-Venugopalan model, which assumes that W^^ [p] is a Gaussian^^'^^'^^ : 



^n, [P] = exp 



/ 



p{x±)p{x±) 



(96) 



The idea behind this model is that the color charge per unit area, p{x±), is the 
sum of the color charges of the partons that sit at approximately the same impact 
parameter. In a large nucleus, this will be the sum of a large number of random 
charges; for Nc = 3, this leads to a Gaussian distribution for p plus a small (albeit 
physically very relevant) contribution from the cubic Casimir The fact that 
this Gaussian has only correlations local in impact parameter is a consequence 
of confinement : color charges separated by more than the nucleon size cannot be 
correlated. The MV model is generally used at a moderately small x, of the order of 
10~^. If the problem under consideration requires smaller values of x, one should use 
the BK or JIMWLK equations, with the MV distribution as the initial condition. 



3.5. Analogies with reaction- diffusion processes 

There are interesting analogies between the evolution equations that govern the 
energy dependence of scattering amplitude in QCD and simple models of reaction- 
diffusion processes^^ . The simplest setting in which these correspondences can be 
seen is to consider the dipole scattering amplitude off a large nucleus, and to assume 
translation and rotation invariance in impact parameter space. It is useful to define 
its Fourier transform as 

r ^ .. (T(0,x±)) 
N{Y,k^) = 2Trj d^xi_ 6'*=^"=^ ^ ^ . (97) 

(Note the factor l/x^ included in this definition.) It turns out that for this object 
N, the BK equation has a very simple non-linear term, 

dN{Y,k^) ^ a,N, 

dY TT 

In this equation, L = ln{kj_/kQ) and x{l) = 2tp{l) — ipij) — — 7) with tp{z) = 
dlnr{z)/dz. The function x{l) has poles at 7 = and 7 = 1, and a minimum 
at 7 = 1/2. By expanding it up to quadratic order around its minimum, and by 
defining new variables, 

t^Y 

z^L+^x"{mY, (99) 
the BK equation simplifies into 

dtN = dlN + N-N'^ , (100) 



xi-dL)NiY,k^)~N\Y,k^) . (98) 
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known as the Fisher-Kolmogorov-Petrov-Piscounov (FKPP) equation. This equa- 
tion has been extensively studied in the literature, because it is the simplest real- 
ization of the so-called reaction-diffusion processes. It describes the evolution of a 
number N of objects that live in one spatial dimension. The diffusion term d'^N 
describes the fact that these entities can hop from one location to neighboring lo- 
cations. The positive linear term +N means that an object can split into two, and 
the negative quadratic term — iV^ that two objects can merge into a single one. One 
can easily check that this equation has two fixed points, = which is unstable 
and iV = 1 which is stable. 

An important property of this equation is that it admits asymptotic travelling 
waves as solutions. Let us assume that the initial condition N{to,z) goes to 1 at 
z —^ — oo and to at z ^ -|-oo, with an exponential tail N{to, z) exp(— /3z). 

z — >+oo 

If the slope of the exponential obeys /3 > 1, the solution at late time depends only 
on a single variable, 

N{t,z) ^ N{z~2t-\\n{t)) . (101) 

When t +00, the logarithm can be neglected in front of the term linear in time, 
and one has a travelling wave moving at a constant velocity dz/dt = 2 without 
deformation (see figure 13). Moreover, this velocity is independent of the details of 




z 



Fig. 13. Travelling wave solutions of the FKPP equation. Red : initial condition. Magenta : solution 
at equally spaced times. 



the initial condition for a large class of initial conditions. 

Going back to the dipole scattering amplitude, this result implies the following 
scaling behavior at large Y : 

{T{0,x^))^^T{Q,iY)x^) , (102) 

with a saturation scale of the form 

Ql{Y) = kl Y^^ e"^ . (103) 

(The exponential comes from the constant in the velocity of the travelling wave, 
and the power law correction comes from the subleading logarithm.) This scaling 
property has an interesting phenomenological consequence for the inclusive DIS 
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Fig. 14. Photon-proton total cross-section measured at HERA, displayed against t = /Qs{Y). 



cross-section, that one can express in terms of the forward dipole scattering ampH- 
tude thanks to the optical theorem : 



^1 p 



(104) 



In this formula, if'i^i ^-Lj Q^) is the light-cone wave function for a photon of virtu- 
ality that splits into a quark-antiquark dipole of size x±, the quark carrying 
the fraction z of the longitudinal momentum of the photon. This wavefunction can 
be calculated in QED, and its only property that we need here is that it depends 
only on the combination [m^ -|- Q^z^{l — z)^]x\ where m is the quark mass. If 
one neglects the quark mass, then eq. (102) implies a simple scaling for the j*p 
cross-section itself : 



o{Y,Q^)=a^,p{QyQl{Y)) 



(105) 



Such a geometrical scaling^^ has been found in the DIS experimental results'^, as 
shown in figure 14. A comment is in order here; as the approach based on coUinear 
factorization and the DGLAP equation succeeds at reproducing much of the inclu- 
sive DIS data, it certainly also reproduces this scaling that is present in the data. 
However, this approach does not provide an explanation for the scaling. It arises 
via some fine tuning of the initial condition for the DGLAP evolution. In contrast, 
in the Color Glass Condensate description of DIS, this scaling is almost automatic. 



"^In addition to explaining geometrical scaling, saturation inspired fits of DIS data are quite suc- 
cessful at small x. See 
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4. Lecture III : Nucleus-nucleus collisions in the CGC framework 
4.1. Introduction 

Up to now, we only considered DIS, in which a possibly saturated proton or nucleus 
is probed by an elementary object'' - a virtual photon that has fluctuated into a 
quark-antiquark dipole. In such a situation, the scattering amplitude can be written 
in closed form as a product of Wilson lines, and its energy dependence can be 
obtained either from Balitsky's equations or from the JIMWLK evolution of the 
distribution of sources that produce the color field of the proton. There are however 
interesting problems that involve two densely occupied projectiles. The archetype 




Fig. 15. Typical contributions to gluon production in hadronic collisions. The dots denote the 
color sources. Left: dilute regime. Right: saturated regime. 

of such a situation is a high-energy nucleus-nucleus collision. In these collisions, 
one of the main challenges is to calculate the multiplicity of the particles (gluons 
at leading order) that are produced at the impact of the two nuclei. In the Color 
Glass Condensate framework, one has to couple the gauge fields to a current that 
receives contributions from the color sources of the two projectiles, 

= S^+S{x-)piix^) + S^~Six+)p2{x^) . (106) 

The fact that there are two strong sources leads to complications that are two-fold : 

• there is no explicit formula that gives the multiplicity (or any other observ- 
able) in terms of Wilson lines in the collision of two saturated projectiles, 

• if one is interested by the particle spectrum at some rapidity Y , one must 
evolve the two projectiles from their respective beam rapidity to Y. The 
question of the factorization of the large logarithms of 1/x is now much 
more complicated than in DIS. 

The kind of complications one is facing in this problem is illustrated in figure 15. In 
the saturated regime, reactions initiated by more than one parton (color source in 
the CGC description) in each projectile become important. Moreover, there can be 



'^Proton-nucleus collisions also belong to this category. Examples of processes have been studied 
in^*. 
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a superposition of many independent scatterings, that will appear as disconnected 
graphs. 



4.2. Power counting and bookkeeping 

In the saturated regime, the color density p (represented by dots in figure 15) is non- 
perturbatively large p ^ g~^. This is due to the fact that the occupation number, 
proportional to (pp), is of order aj^ in this regime. Thus for a connected graph, 
the order in g is given by 



- g^^r. , (107) 



where Ug is the number of produced gluons and the number of loops. One can see 
that this formula is independent of the number of sources p attached to the graph. 
Indeed, since each source brings a factor g~^ and is attached at a vertex that brings 
a factor g, each source counts as a factor 1. If the diagram under consideration is 
made of several disconnected subgraphs, one should apply eq. (107) to each of them 
separately. 

Among all the diagrams that appear in the calculation of particle production, 
a special role is played by the so-called vacuum diagrams - diagrams that have 

rig = external gluons. They only connect sources of the two projectiles, and 
are thus contributions to the vacuum-to- vacuum amplitude (OoutlOin), hence their 
name. The order of connected vacuum diagrams is gf^^"^"^^. An extremely useful 
property is that the sum of all the vacuum diagrams (connected or not) is the 
exponential of those that are connected (that we denote iV[j] where j is the external 
current due to the color sources of the two projectiles) 

E /all the vacuum \ / connected \1 ,T/r,i 

= exp O )} =e" . (108) 

\ diagrams / V vacuum diagrams/ J 

The reason why vacuum diagrams are important in our problem is that it is possible 

to write all the time ordered products of fields - that enter in the reduction formulas 
for gluon production amplitudes - as derivatives of exp{iV[j]) 

(0out|TA(xi) • • • A(x„)|Oi„) = — ^ • • • — ^ e^^IJ'l . (109) 
' ' ^SJ{xl) iS j{xn) 

Thanks to this property, one can write a very compact formula for the probability 
Pn of producing exactly n gluons in the coUision^^ 



. . . > (110) 

J+=3-=J 
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where the operator T> is defined by* 

■r./^o»_(.,riaD.j^^ , ^^^^^ 

An important point to keep in mind about cq. (110) is that the external currents 
must be kept distinct in the amplitude and complex conjugate amplitude until all 
the derivatives contained in T) have been taken. Only then one is allowed to set 
j+ and j_ to the physical value of the external current. The propagator G\_, 
that has only on-shell momentum modes, is the usual cut propagator that appears 
in Cutkosky's cutting rules^^'^^. The operator V acts on cut vacuum graphs by 
removing two sources (one on each side of the cut, i.e. a j+ and a j-), and by 
connecting the points where they were attached by the cut propagator G'^_. In fact, 
since P„ is obtained by acting n times with the operator V, it is the sum of all the 
cut vacuum diagrams in which exactly n propagators are cut. Eq. (110) also makes 
obvious the fact that the probabilities P„ do not have a meaningful perturbative 
expansion in the saturated regime, because the sum iV\j] of the connected vacuum 
diagrams starts at the order g~'^. 

By summing eq. (110) from n = to oo while keeping j+ and j- distinct, 
one obtains the sum of all the cut vacuum diagrams with the current j+ in the 
amplitude and j- in the complex conjugate amplitude to be 

y/ all the cut \ ^ ^^^^^ 
^-^ \ vacuum diagrams / 

When we set j+ = j- , this sum becomes J2n therefore it should be equal to 

1 because of unitarity. Eq. (110) is very useful, because it allows to replace infinite 
sets of Feynman diagrams by simple algebraic equations. Similarly, the fact that 
eq. (112) is 1 when j+ = j_ corresponds to a cancellation of an infinite set of 
graphs", that would be very difficult to see at the level of diagrams. 



4.3. Inclusive gluon spectrum 

Eq. (110) leads to compact formulas for moments of the distribution of produced 
particles. The first moment - the average multiplicity - reads^^ 

C30 

N = V n P„ = P I e^^[^+l e"*^*!^-!) . (113) 

With the help of eq. (112), this formula tells us that N is given by the action of 
the operator V on the sum of all the cut vacuum diagrams. In plain english, this 

*We are a bit careless here with the Lorentz indices, polarization vectors, etc, because our main 
goal is to highlight the general techniques for keeping track of the diagrams that contribute to 
particle production in the saturated regime. 

"This cancellation is closely related to the AbramovsJcy-Gribov-KancheJi canceJJation^^ . 
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translates into : take a cut vacuum diagram (connected or not), remove a source 
on each side of the cut, and put a cut propagator where the sources were attached. 
Depending on whether the cut vacuum diagram one starts from is connected or not, 
one gets two different topologies, displayed in figure 16. Each of the blobs in these 




Fig. 16. The two topologies contributing to the average gluon multiplicity N. In each blob, one 
must sum over all the possible ways of cutting the propagators. 




diagrams can be any connected graph, and must be cut in all the possible ways^. 
Thus, only connected graphs contribute to the multiplicity. 

An important point is that, even though the perturbative expansion for the 
P„ is not well defined, the multiplicity (and more generally any moment of the 
distribution P„) can be organized in a sensible perturbative series^. The Leading 
Order is obtained by keeping only the leading order vacuum graphs, i.e. those that 
have no loops : ^ 

^.o = E E dj I i';^- (114) 

trees cuts 



Thus N starts at the order g^^. In eq. (114), for each tree diagram, one must sum 
over all the possible ways of cutting its lines. The simplest way of doing this is to 
use Cutkosky's rules : 

• assign + or — labels to each vertex and source of the graph, in all the 
possible ways (there are 2" terms for a graphs with n vertices and sources) . 
A + vertex has a coupling —ig and a — vertex has a coupling +ig, 

• the propagators depend on which type of labels they connect. In momentum 
space, they read : 

G'^^{p) = i/ip^ + if) (standard Feynman propagator) 

G^--{p) = ^^l{v^ ^ (complex conjugate of G° _|_(p)) 

G"+„(p) = 2^(?(V)'5(p') 

G"_+(p) = 2^(?(pO)<5(p2) . (115) 

^Note that by not performing the d^p integration contained in the explicit cut propagator, one 
obtains the inclusive gluon spectrum dN /cP'p instead of the integrated multiplicity. 
"'The fact that this is possible for A'' but not for the Pn's themselves is due to the fact that the 
only graphs that contribute to N are connected. This is a consequence of the AGK cancellation. 
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A quick analysis shows that, when one sets j+ — j_, summing over the ± labels at 
each vertex produces combinations of propagators, 

G%(p)-G°_(p) = G°(p) 

G"_+(p)-G°__(p) = G°(rt, (116) 

where G° {p) is the retarded propagator^. Thus, for a given tree graph, doing the 
sum over the cuts simply amounts to replacing all its propagators by retarded 
propagators. The last step is to perform the sum over all the trees. It is a well 
known result that the sum of all the tree diagrams that end at a point a; is a 
solution of the classical equations of motion of the field theory under consideration. 
In our case, this sum is a color field A^{x) that obeys the Yang- Mills equations 

[2?^,^^ = J^^ (117) 

where J" is the color current associated to the sources pi,2 that represent the 
incoming projectiles (see eq. (106)). The boundary conditions obeyed by A'^{x) 
depend on the nature of the propagators that entered in the sum of tree diagrams. 
When these propagators are all retarded, one gets a retarded solution of the Yang- 
Mills equations, that vanishes in the remote past, limj,o_>_oo A^{x) = 0. The precise 
formula for the gluon spectrum in terms of this solution of the Yang- Mills equations 
reads 



^ j d^xd^y e^f(--^) U,Uy Y,4^l M^)My) ■ (118) 



dYd'^p^ 167r3 

Note that, although the integrations over x and y look 4- dimensional, they can be 
rewritten as 3-dimensional integrals evaluated at xq — > +oo, thanks to the identity 

/ d'^x e'f^ Af,{x) = lim / d^x e'f=" [^o - lEp] ^^(x) . (119) 

J x°^+oo J 

Solving the Yang-Mills equations is an easy problem in the case of a single soiircc p, 
but turns out to be very challenging when there are two sources moving in opposite 
directions. The Schwinger gauge, defined by the constraint A'^ = x~^A~ +x~A~ = 
0, is quite useful because it alleviates the need to ensure that the current J'' is 
covariantly conserved''. In this gauge, A'^ ~ where J~ ^ and conversely, which 
makes this condition trivial. Moreover, in this gauge, one can find the value of 
the gauge field on a time-like surface just above the light-cone (at a proper time 
r = 0+) simply by matching the singularities across the light-cone. These initial 



^In momentum space, G'^(p) = i/{p^ + isign(po) e)- Therefore, in coordinate space, it is propor- 
tional to ^(a;" — y"), hence its name. 

''In general gauges, one has to enforce the condition [O^, J''] = (this is a consequence of 
Jacobi's identity for commutators). Because this relation involves a covariant derivative rather 
than an ordinary derivative, the radiated field leads to modifications of the current itself. 
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conditions'^ can be written as^ 



A^{t = 0,x^) = -| [A\{x^) , AUx±)] 



= (gauge condition) , 



(120) 



where A^ = t ^(.t A'^ — x A ). In this formula, Ai{xx) and A2{x±) are the 
gauge fields created by each nucleus'^ below the light-cone : 



Therefore, the problem of solving the Yang-Mills equations from xq = — oo to 
xq = +O0 is reduced to solving them in the forward light-cone from a known initial 
condition''. 

Since our problem is invariant under boosts in the z direction, one can com- 
pletely eliminate the space-time rapidity t] from the equations of motion (and the 
initial conditions in eq. (120) are also ry-independent). Thus, in the forward light- 
cone, one has to solve numerically'^ equations of motion in time and two spatial 
dimensions, and then to evaluate cq. (118). The result of this computation is dis- 
played in figure 17. In this computation, the MV model was used as the distribution 
of the sources pi and p2. Therefore, the dependence of the spectrum on the mo- 
mentum rapidity Y of the produced gluon cannot be obtained in this calculation, 
and only the k± dependence is shown. The main effect of gluon recombinations on 
this spectrum is that it reduces the yield at low transverse momentum, k± ^ Qs- 
Indeed, in a fixed order calculation in perturbative QCD, the spectrum would be- 
have as fcj^. In the CGC picture, the singularity of the spectrum at low k± is only 
logarithmic'^, and is therefore integrable. 

4.4. Inclusive quark spectrum 

A similar study has also been performed for the initial production of quarks in 
nucleus-nucleus collisions'''. The starting point is to construct for quarks an oper- 

''An interesting feature of the gauge fields at early times after the collision — a phase recently 
named "glasma" - is that the chromo-electric and magnetic fields are purely longitudinal, while 
they were transverse to the beam axis just before the coUision^^. 

''Because retarded solutions are causal, the field below the light-cone cannot depend simultaneously 
on pi and p2- 

''Note that at r > 0, the YM equations are the vacuum ones, since all the sources are located on 
the light-cone. 

"^If the final Fourier decomposition is performed at a finite time r, the spectrum is completely 
regular when k± —* 0. 



A 



A\ 



■{ = -Ui{x^)d'ul{x^) 



i = -U2{x^)d'uUx^) 
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Fig. 17. The gluon spectrum at leading order in the CGC framework. 



ator Vq that plays the same role as the operator V defined in eq. (Ill) : 

T^,^ f Sl_{x,y)0Jy \ J (122) 
Jx,y Sri+{x) 6ri_{y) 

where S^_{x, y) is the free cut fermionic propagator and where 77 is a Grassmanian 
current that couples to the spinors. In terms of this operator, the probability of 
producing n quarks is given by : 



.... (123) 

The first thing to note is that now the connected vacuum diagrams, whose sum is 
iV ^ depend on both the source j and on the source 77. However, the latter is set to 
zero at the end of the calculation, because in the CGC one assumes that the color 
sources in the wavefunction of the projectiles couple only to the gluons. Therefore, 
the source r\ serves only as an intermediate bookkeeping device. Another important 
point in this formula is the presence of the factor exp(I?). This factor means that 
we are calculating an inclusive probability, for producing exactly n quarks possibly 
accompanied by an arbitrary number of gluons*^. In practice, this fact means that 
one must sum over all the possible ways of cutting the gluons lines in the diagrams 
that contribute to quark production. From eq. (123), one obtains the following 
formula for the average number of produced quarks 



■ (124) 

J4-=J_=J ^ ^ 

T7^ = T7_=0 



In this formula, the underlined factors represent the sum of all (connected or not) 
the cut vacuum diagrams made of quarks and gluons, with sources ^+,774- on one 



"^Without this factor, we would be calculating the probabiUty of producing n quarks and gluons. 
Note that in principle, we should also modify our definition of the probability of producing n 
gluons by a factor exp(I'q). However, the quarks are a subleading correction compared to the 
gluons, and this change would not affect the gluon spectrum at leading order. 
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side of the cut, and sources j-,r]- on the other side. Acting on a term of this 
sum with Vq removes a source 77+ and a source r^-, and connect the points where 
these sources were attached by a cut fermion propagator. Diagrammatically, this 
corresponds to the two topologies displayed in figure 18. Note however that the 




Fig. 18. Topologies corresponding to oq. (124). 



topology that appears on the left of figure 18 cannot exist because it has a quark 
line which is not closed onto itself (this is forbidden since we set the fermionic 
sources r]± to zero at the end of this calculation). Thus, we only have the second 
family of diagrams, that have at least one loop. This means that the average number 
of quarks is of order g", compared to the number of gluons which is of order g^^. 

The leading contribution to the quark multiplicity is obtained by including only 
tree diagrams in the blob. Thus, we have to sum all the graphs that have one quark 
loop (with an explicit cut on it) and any number of gluonic trees attached to it, 
and all the cuts thereof. The sum of all the gluonic trees and their cuts has already 
been encountered in the computation of the gluon multiplicity : it is equal to the 
retarded solution A'^ix) of the Yang-Mills equations that vanish in the remote past. 
Therefore, the quark spectrum is given by 



dNq _ 1 



dYdPp^ 167r3 



e^P-^u{p) S+-{^,y) 0^ u{p)e-'P-y , (125) 



y 



where 5-| is the cut quark propagator on which the retarded classical field A'^ has 

been resummed. This resummed propagator can be obtained as the solution of the 
equation 

S,A^,y)^Si,{x,y)~igJ2i-'^r f d'zS%ix,z)A^iz)rSrjAz,y) , (126) 

where e,e' = ± (we need only the combination e = -|-,e' = — in eq. (125), but 
the four terms get mixed when one resums the background field). It is possible to 
decouple these equations by performing a "rotation" on the e, e' indices^^. 

See' — * Sa(3 = ^ UaeUpe'See' 

i-iySee' ^ = ^ {/„,{/^,(-l)' , (127) 



e=± 
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With f/=;^ (;-;)• 

After this rotation, the propagator matrix becomes triangular, 

with Sj^ and the resummed retarded and advanced propagators and where 
S^{p) = 2tt-^S{p'^). The main simpHfication comes from the fact that the product of 
the free matrix propagator and of S is the sum of a diagonal and a nilpotent matrix, 
which makes the calculation of its n-th power very easy**. In particular, one finds 
that the equations that lead to the retarded (and also the advanced) propagator do 
not mix with anything else, 

S^{x,y) = Sl{x,y)-ig j d'z Sl{x,z)A^{z)^^S^{z,y) , (131) 

and that the resummed 5^ can be expressed in terms of as^ 

5, =5,* 50-1* 50*50-1* 5,. (132) 

At this point, one must invert the rotation done in eq. (127) in order to obtain 

S-i which is needed in the formula for the quark spectrum. This gives the quark 

spectrum in terms of retarded quantities, 

dJVq 1 f d^o 2 



1 f d^q 



where is the "scattering part" of the retarded propagator, related to by 

Sn=Sl + Sl*T,*Sl. (134) 

The last step is to write this object in terms of retarded solutions of the Dirac 
equation in the background field A'^. It is easy to check that 

(p, q) = hm / d^x e^f- (p) V, {x) 

x"— ►+00 J 

{i0,-gMx))xl,^{x) = Q, V,(xO,x) = v{q)e"'-- . (135) 

— *— 00 

In this formula, u{p) and v{q) are the usual free spinors. Note that the initial 
condition for the Dirac equation is a negative energy spinor, and that the projection 
performed at the final time is with a positive energy spinor. In the vacuum, there 
would be no overlap between these spinors. However, since in our problem the 
spinor travels on top of a time-dependent background field, it acquires positive 
energy modes which make non zero. 

''Indeed, with very formal notations, the resummed matrix propagator is 

00 

S = sO^(-i5^r[sS°]". 

n=0 

'The * symbol denotes the convolution of 2-point functions: {A * B){x, y) = A{x, z)B{z, y). 
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Fig. 19. Numerical results on quark production from the CGC. Left: time evolution of the quark 
yield. Right: quark fcx spectra for different masses. 



This formulation of quark production in the CGC framework has been imple- 
mented numerically, also with the MV model for the average over the configurations 
of the color sources pi^2 - Similarly to what happened with the gluons, one can obtain 
analytically the value of the spinors just above the light-cone. Hence, the numerical 
resolution of the Dirac equation is only needed in the forward light-cone. However, 
there is a major difference compared to the gluons at LO : even though the back- 
ground color field does not depend on rapidity, this is not true of the solutions 
of Dirac equation^. Indeed, the momentum q in the initial condition renders the 
spinors dependent on the space-time rapidity 77 (the boost invariance of the back- 
ground field implies that the spinors depend only on the difference rj — yq where 
Uq is the rapidity of the momentum q). This difference makes the computation of 
the quark spectrum much more computationally intensive relative to that of the 
gluon spectrum, because one has to keep the three dimensions of space. Some of 
the results obtained are displayed in figure 19. On the left plot is shown the time 
dependence of the quark yield, for different quark masses (i.e. the yield obtained 
by performing the projection in eq. (135) at a finite time instead of taking the 
limit xq -|-cxd). One can see that a good fraction of the quarks are produced at 
T — 0, when the two nuclei pass through each other** and that the number slightly 
increases in time afterwards due to the color field present in the forward light-cone. 
The right panel of figure 19 shows the k± dependence of the spectrum for various 
quark masses. As expected, the spectrum is harder for larger quark masses. Note 
that the tail of the curves is probably affected by important lattice artifacts due to 
a too coarse lattice. 



^This has nothing to do with the fact that we are considering fermions, but rather with the quark 
spectrum being a NLO quantity - that involves a loop in the background of the classical field. 
''In the analogous QED problem of e+e" production in the high-energy collision of two electrical 
charges, all the electrons are produced at r = and their number does not change at r > 0. This 
is because in QED, the electro-magnetic potential in the forward light-cone is a pure gauge, that 
could be made to vanish by a gauge transformation. 
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4.5. Loop corrections to the gluon spectrum 

Thus far, wc limited ourselves to the leading order contribution for both the gluons 
and the quarks. However, we a priori know from figures 16 and 18 what diagrams 
contribute to the gluon and quark multiplicities to all orders. There is therefore 
a well defined and systematic procedure to compute corrections to the previous 
results. Loop corrections to gluon production are very relevant for the following 
reasons: 

• They contain terms that are divergent due to unbounded integrals over 
longitudinal momenta, very similar to the divergences encountered in the 
derivation of the BK equation. One should verify whether these divergences 
can be absorbed in the distributions W^[pi] and VF[p2] of the color sources 
of each projectile. This factorization is crucial for the internal consistency 
of the CGC framework. 

• It has been noted recently that the boost invariant solution A'^ix) of the 
Yang-Mills equations is unstable'; rapidity dependent perturbations to this 
solution grow exponentially in time. Loop corrections generate this kind of 
rapidity dependent perturbations. Tracking all these terms and resumming 
them is very important in order to get meaningful answers from the CGC 
regarding the momentum distribution of the produced gluons, and may be 
relevant in the problem of thermalization in heavy ion collisions. 

Note that these two items address very different stages of the collision process. The 
first relates to the incoming wavefunctions (and as such should be independent of 
the subsequent collision), while the second issue is about what happens in the final 
state after the collision. Therefore, we should aim at writing the 1-loop corrections 
in a way that separates the initial and final state as clearly as possible. 

Let us start by listing the relevant diagrams : the 1-loop corrections to the 
average multiplicity are shown in the diagrams of figure 20. The topology on the 




Fig. 20. 1-loop diagrams contributing to the gluon spectrum. 



left is very similar to the one already encountered at tree level, except that one of 
the blobs has now a loop correction in it. The topology on the right is new; but it is 
in fact similar to what we had to evaluate in the case of the quark spectrum, except 



'This instability is very similar to the Weibel instability that occurs in anisotropic plasmas' 
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that the fermionic cut propagator must be replaced by the cut propagator 

G+- of a gluon. The NLO contribution to the gluon spectrum can be written as 

167r3 



dYdPp^ 



— I 

67r3 J 



{x)5A' {y)+5A^' (a;)^" (y) 



+Gr-(.x,y) 



(136) 



The two terms of the first line are the contribution of the diagram on the left of 
figure 20 (the loop can be in either of the two blobs), and the term on the second 
Unc comes from the diagram on the right. The field 5 A that appears on the first 
line is the 1-loop correction to A; and it obeys the linearized equation of motion 
for small fluctuations. 

Let us now illustrate how one can separate the initial state from the final state 
in the term that contains G'i^_ {x, y). First, by analogy with the case of the quarks, 
we can write 

d^q 



I 



d xd ye 



ip-{a 



A A, A' 



lim 



d^x e^f^ {dl 



X'q 



(137) 



where a^/q(a;) is a small fluctuation of the gauge field on top of A'^, with initial 



condition e't,e^'^'^ when xq 



-00. The equation of motion of this fluctuation is 



obtained by writing the Yang- Mills equations for ^ + a and by linearizing it in a. 
A central formula in order to separate the initial and final states is the following-i 



a{x) = j d^y 

r=0+ 



a{Q,y)-Ty A{x) 



(138) 



where (0, y) denotes a point located on the light-cone (t = 0) (y represents any set 
of three coordinates that map the light-cone.) In this formula, the classical field A 
is considered as a functional of its initial condition ^(0, y) on the light-cone. The 



notation 



a{Q,y)-Ty 



is a shorthand for 



a{Q,y)-Ty =a{0,y) 



SA{0,y) 



+ 



(n • dy)a{0,y) 



S{n-dy)A{0,y) 



(139) 



(In this formula, the 4- vector n'' is a vector normal*^ to the light-cone.) The proof of 
eq. (138) is straightforward', but its diagrammatic interpretation is more interesting. 
Note first that ^''(a;), seen as a functional of its initial condition on the light-cone, 
can also be represented by tree diagrams, as illustrated in the left panel of figure 21. 
(This can be seen from the Green's formula for A{x).) The action of the operator 



J To avoid encumbering the equations with indices of various kinds, we ajre suppressing all the 
indices in this and the following formula. 

^ni^dx^ = if dXfx is a small displacement on the light-cone at the point under consideration. 
'Write the Green's formula that expresses A{x) in terms of the initial A{0, y), insert it in eq. (138), 
and check that this leads to the Green's formula that relates a{x) to its initial condition a{0,y). 
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AW 





Fig. 21. Left : diagrammatic representation of ^ as a function of its initial condition on the light- 
cone (the open dots represent the initial .4(0, y)). Right : propagation of a small fluctuation on 
top of the classical field. 



defined in eq. (139) on the classical field A{x) is to replace one of the open dots in 
figure 21 by the fluctuation a{0,y), represented by a filled dot in the right panel 
of figure 21. The diagram one gets after this is nothing but a contribution to the 
propagation of a small fluctuation over the classical field. Plugging eq. (138) in 
eq. (137), this quantity becomes 



lim 



X 
r=0+ 



d^xd^y e'P-^^-y^ (9° 
d\ 



zEp){d^ + zEp)J24^t 



aA'q(0,M) -T^] ^"(x)] [[a^,q(0,-i;) -T^] A''{y) 



(140) 



The brackets are crucial in this formula, in order to limit the scope of the derivatives 
contained in the operators and T^. Note that, if it were not for these brackets, 
the first line and the two ^'s of the second line would be nothing but the LO 
gluon spectrum. It turns out that, after one adds the contribution of the first line 
in eq. (136), the NLO correction to the spectrum can be written as 



dN. 



dYd^p^ 



5A{Q, u) ■ 



d-^ud V 



.T = 0+ 



T = + 



dN, 



dYd^Pj_ ' 
(141) 

where the LO spectrum is considered as a functional of the initial classical field on 
the light-cone. In this equation, 6A{0, u) is the value of SA on the light-cone, and 
the 2-point S(m, v) is defined as 

d^q 



aA'q(0,u)a^,q(0,'u) . 



(142) 



Note that M(0, m) and Yj{u,v) arc in principle calculable analytically. Eq. (141) 
realizes the separation we were seeking of the initial and final states. Indeed, the 
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Fig. 22. Illustration of eq. (141). The 1- and 2-point functions below the light-cone are respectively 
Sji{0, u) and v). 

operator in the square bracket depends only on what happens below the light-cone, 
i.e. before the collision. On the contrary, the LO spectrum seen as a functional of 
the initial classical field A depends only on the final state dynamics. The other 
benefit of this formula is that is expresses the NLO correction as a perturbation of 
the LO one; this property ~ that seems generalizable to other inclusive observables 
- suggests the universality of the initial state divergences and their factorizability. 

From eq. (141), it is easy to see what are the potential sources of divergences. 
A first issue is that the coefficients 6A{0,u) and v) are infinite. For 
for instance, the integration over the longitudinal component of the momentum q 
in eq. (142) diverges. A similar divergence occurs in the loop contained in SA{0, u). 
The fact that these divergences arise in the first factor of eq. (141) indicates that 
they are related to the evolution of the initial projectiles prior to the collision. 
These divergences can be momentarily regularized by introducing cutoffs Yq , Yq in 
rapidity around the rapidity Y at which the spectrum is calculated. Thus, SA{0, u) 
and E(u, d) become finite but depend on these unphysical cutoffs. To be consistent, 
the distribution of the sources pi and p2 should be evolved from the beam rapidities 
to Yq and Yq respectively. Thus, the complete formula for the LO-f NLO spectrum, 
including the average over the sources, should be 



dYd^p 



LO+NLO 



I [Dp,][Dp,] W^.,_-.„[Pi]W^.,_,.JP2] 



1 ^0 




u SA{0, u)-Tu + / d^ud^v E(m, v) ■ T^T.. 




V 



dYd^p^ ' 




(143) 



where the subscript yj and superscript Yq indicate that the momentum integrals 
contained in the bracket have cutoffs in rapidity. Recall that the LO spectrum in 
the right hand side is a function of A on the light-cone, which is itself a function 
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of pi,2- The factorizability of these divergences in the initial state is equivalent to 
the independence of the previous formula with respect to the unphysical Yq and Yq. 
Let us for instance change Yq into Yq + AYq. According to the JIMWLK equation, 
the distribution of pi is modified by 

^v-,_-v-o[Pi] - [l + AYon[pi]]w,^__,^[pi]. (144) 

At the same time, the operator in the right hand side varies by 

0^°,[pi,P2] - Ol°,[pi,p2]+AYo^^ . (145) 

At this point, one can verify that the terms hnear"^ in AYq cancel provided that 

dO^°,[puP2] 

=Ht[pi]. (146) 

Similar considerations on the Yq dependence give another condition : 



dOl",[puP2] 



dY^ 



= -n^[P2]. (147) 



Therefore, in order to check whether one can factorize these divergences in the 
JIMWLK evolution of the incoming sources, one must calculate the coefficients 
dA{0,u) and T,{u,v) - or at least their divergent part - and remap the operator 
0^,[pi,P2] into the JIMWLK Hamiltonian. Although this program has not been 

fully implemented yet, one can already note that the structure of O^, [pi, /02] makes 
this outcome very plausible. 

Eq. (141) also allows us to discuss the issue of the instability of the boost 
invariant classical solution. These instabilities manifest themselves in the fact that 
the action of r„ on A{x) diverges when the time xq goes to infinity. Indeed, 

is a measure of how A{x) is sensitive to its initial condition. Therefore, if the solution 
A{x) is unstable, small perturbations of its initial condition lead to exponentially 
growing changes in the solution. Prom the numerical study of these instabilities (see 
figure 23), one gets^^ 

T^A{x) ~ , (149) 

where p is of the order of the saturation momentum. This means that, although 
the 1-loop corrections are suppressed by a factor as compared to the LO, some 
of these corrections are enhanced by factors that grow exponentially in time after 



™Note that, since we have only considered 1-loop corrections, this independence can only be 
satisfied for small variations of the cutofi', at lineaj: order in these variations. 
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Fig. 23. Time dependence of small fluctuations on top of the boost independent classical field. 



the collision. At first sight, one may expect a complete breakdown of the CGC 
description at 



W-Q^Mn^ i^-Lj ^ (150) 

i.e. the time at which the 1-loop corrections become as large as the LO contribution. 
The only way out of this conclusion is to resum all these enhanced corrections in 
the hope that the resummed series is better behaved when r —^ +00. Let us assume 
for the time being that we have performed this resummation, and that the sum of 
these enhanced terms generalize eq. (141) to read 

dY<Pp^ ^^^'^^ dY<Pp^ ' ^^^^^ 

where Z\Tu] is a certain functional of the operator T^. In the right hand side, 
we have emphasized the fact that the LO spectrum is a functional of the initial 
classical field on the light-cone. This formula can be written in a more intuitive 
way by performing a Fourier transform of Z[T„], 

Z[T„] = j [Da{u)\ e*^^=o+^'" [aW-T„] z[a{u)] . (152) 

In this formula, the functional integration [Da{u)\ is in fact an integration over 
two fields : the fluctuation a{u) itself and its derivative normal to the light-cone 
{n ■ du)a{u). Thanks to the fact that Tu is the generator of translations of the 
initial conditions on the light-cone, the exponential in the previous formula is the 
translation operator itself. Therefore, when this exponential acts on a functional 
of the initial classical field yt(0, m), it gives the same functional evaluated with a 
shifted initial condition A{Q,u) + a{u). Therefore, we can write 

dlV_ed ^ f ^^^^^^^ -j^^^^j dN ,,[A{O u) + a{u)] ^^^3^ 



dYd^p^ J ^ ' 'J ' ^ " dVd^p^ 
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We see that the effect of the resummation is simply to add fluctuations to the initial 
conditions of the classical field, with a distribution that depends on the details of 
the resummation". It is easy to understand why these fluctuations are crucial in the 
presence of instabilities : despite the fact that they are suppressed by an extra power 
of as, the instabilities make them grow and eventually become as large as the LO. 
One can also see that the resummation has the effect of lifting the time limitation 
of cq. (150). Indeed, after the resummation, the fluctuation a{u) has entered in the 
initial condition for the full Yang-Mills equation, whose non-linearities prevent the 
solution from blowing up. A very important question is whether these instabilities 
fasten the local thermalization of the system formed in heavy ion collisions. 



4.6. Summary and outlook 

If the initial state factorization works as expected, and after the resummation of the 
leading contributions of the instability, the formula for the gluon spectrum should 
read 



dN 



/[ 



Da] m . (154) 

This formula resums the most singular terms at each order in ag ■ Because of their 
relation to the physics of the initial and final state respectively, the distributions 
W[p] generalize parton distributions, while Z[a] plays a role similar to that of a 
fragmentation function". 

Note that, even after the resummations performed in the initial and final states 
of eq. (154), this formula still suffers from the usual problem of collinear gluon 
splitting in the final state. This is not a serious concern in heavy ion collisions 
though, because collinear singularities occur only when one takes the r +oo limit, 
and we do expect to have to switch to another description (like hydrodynamics) 
long before this becomes a problem. In fact, the initial condition for hydrodynamics 
should be specified in terms of the energy- momentum tensor, which is infrared and 
collinear safe because it measures only the flow of energy and momentum. 

A more important problem, that has still not received a satisfactory answer, is 
to understand how the initial particle spectrum or the local energy momentum- 
tensor - become isotropic. This requires formulating a kinetic theory of the glasma 
which describes how particles emerge from the decaying classical field and their 
subsequent interactions both with the classical field and with other particles. Re- 
cently, such a kinetic equation has been derived for a scalar field theory coupled to 



"In a recent work by one of the authors, using a completely diflferent approach, the spectrum of 
initial fluctuations was found to be Gaussian*^ . 

"Naturally, this function has nothing to do with a gluon fragmenting into a hadron. Instead, it is 
related to how classical fields become gluons. 
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strong sources . Extending this work to QCD and exploring its consequences- in 
particular, the approach of the particle+field system towards equilibration remains 
a challenging problem. 
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